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Inhomogeneous "longitudinal" plane waves 
in a deformed elastic material. 
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00 ■ By definition, a homogeneous isotropic compressible Hadamard 

material has the property that an infinitesimal longitudinal homoge- 
q— , neous plane wave may propagate in every direction when the material 

O ' is maintained in a state of arbitrary finite static homogeneous defor- 

• ■ mation. Here, as regards the wave, 'homogeneous' means that the 

Cd i direction of propagation of the wave is parallel to the direction of 

^ ' eventual attenuation; and 'longitudinal' means that the wave is lin- 

early polarized in a direction parallel to the direction of propagation. 
P5 , In other words, the displacement is of the form u = n cos A;(n • x — ct), 

O ' where n is a real vector. 

It is seen that the Hadamard material is the most general one for 

which a 'longitudinal' inhomogeneous plane wave may also propagate 

h>. I in any direction of a predeformed body. Here, 'inhomogeneous' means 

^D ' that the wave is attenuated, in a direction distinct from the direction 

of propagation; and "longitudinal" means that the wave is elliptically 
^Sj I polarized in the plane containing these two directions, and that the 

^ ' ellipse of polarization is similar and similarly situated to the ellipse 

C^ ' for which the real and imaginary parts of the complex wave vector 

are conjugate semi-diameters. In other words, the displacement is of 
the form u = 5R{S exp iuj{S ■ x — ct)}, where S is a complex vector (or 
bivector). 

/\ . Then a Generalized Hadamard material is introduced. It is the 

^ I most general homogeneous isotropic compressible material which al- 

lows the propagation of infinitesimal "longitudinal" inhomogeneous 
plane circularly polarized waves for all choices of the isotropic direc- 
tional bivector. 

Finally, the most general forms of response functions are found for 
homogeneously deformed isotropic elastic materials in which "longi- 
tudinal" inhomogeneous plane waves may propagate with a circular 
polarization in each of the two planes of central circular section of 



the B^-ellipsoid, where B is the left Cauchy-Green strain tensor cor- 
responding to the primary pure homogeneous deformation. 

1 Introduction. 

Infinitesimal longitudinal homogeneous plane waves play a special role in 
classical linearized theory. For such waves the amplitude is parallel to the di- 
rection of propagation n so that all the particles oscillate along the direction 
n. Such waves may propagate in every direction in an isotropic compress- 
ible elastic body. This is no longer the case for propagation in an elastic 
anisotropic crystal. Possible directions of propagation of longitudinal homo- 
geneous plane waves, called 'specific directions' by Borgnis pj, may be as 
few as three in an elastic anisotropic crystal. Of course, by assuming certain 
restrictions on the elastic constants, it is possible, as Hadamard [3] did, to 
have a special model anisotropic material in which longitudinal waves are 
possible in every direction. 

A similar situation holds for infinitesimal motions of an isotropic homo- 
geneous elastic material held in an (arbitrary) state of static (finite) homo- 
geneous deformation - specific directions are limited in number. However, 
it is possible to place conditions on the basic strain energy density, or on 
the response coefficients defining the material, such that every direction is 
specific for arbitrary (finite) homogeneous deformation. The resulting model 
was called a 'Hadamard material' by John [2]. 

Here we consider the propagation of infinitesimal inhomogeneous plane 
waves in an isotropic homogeneous compressible elastic body held in an ar- 
bitrary state of finite static homogeneous deformation. Such waves may be 
described in terms of bivectors - complex vectors - the amplitude bivector A 
and the slowness bivector S, which is written S = A^C, where the directional 
bivector C is written C = mm + ih (m ■ fi = 0, ?7i > 1, |m| = |n| = 1). Once 
the directional bivector C is prescribed, the slowness S and the amplitude 
A are determined from the equations of motion. Prescribing C is equivalent 
to prescribing an ellipse with principal semi-axes mm and n; this direc- 
tional ellipse for inhomogeneous plane waves corresponds to the direction of 
propagation n for homogeneous plane waves [1]. We borrow the adjective 
'longitudinal' to describe an inhomogeneous plane wave for which A and S 
(and C) are parallel: A x S = 0. What this means is that [S] the ellipses 
of A and of S (and of C) are all parallel, similar (same aspect ratio) and 
similarly situated (parallel major axes and parallel minor axes). We deter- 
mine the most general form of strain-energy density so that such waves may 
propagate for all choices of the directional bivector C. It turns out to be a 



Hadamard material. 

Then we consider the possibihty of having 'circularly polarized longitu- 
dinal' waves. For such waves both C and A are isotropic and coplanar: the 
corresponding ellipses are coplanar circles. Such waves are a subclass of lon- 
gitudinal waves. The corresponding class of materials for which circularly 
polarized longitudinal waves are possible is called a 'Generalized Hadamard 
material'. Clearly they include Hadamard materials. 

Next we obtain constitutive equations for materials which allow the prop- 
agation of two special infinitesimal longitudinal circularly polarized inhomo- 
geneous plane waves for all choices of the basic static homogeneous deforma- 
tion. The circles of polarization are to lie in the planes of central circular 
section of the B^-ellipsoid: x-B"x = 1, where B is the left Cauchy-Green 
strain tensor corresponding to the basic deformation. In the particular cases 
where n = ±1, we determine the corresponding forms of the strain-energy 
function. In every case, the slownesses are presented explicitly. Finally we 
present the very simple results when n = |. 

2 Inhomogeneous longitudinal plane waves. 

Here we recall some basic properties of inhomogeneous plane waves and then 
introduce inhomogeneous longitudinal plane waves. Finally we introduce 
inhomogeneous longitudinal circularly polarized plane waves. 

The displacement field corresponding to an infinite train of infinitesimal 
inhomogeneous plane waves is eu (where e is an infinitesimal constant: e^ ^ 
|e|) and u is given by 

u = [A exp 2a;(S ■ X — t)]^ 
= {A+cosu;(S+-x-t) - A"sinw(S+-x-t)}exp-wS" ■ x. (2.1) 

Here A = A"*" + iA^ is the amplitude bivector; S = S^ + 28^ is the slowness 
bivector; and the period of the waves is 27i/u. The particle initially at x is 
displaced to x + eu at time t, so that it moves on an ellipse, centre x, with 
conjugate radii eA"^ exp —ujS~ ■ x and eA~ exp —uS^ ■ x. When the particle 
is at X + eA+exp —uS^ ■ x it is moving parallel to A~, and when it is at 
X -|- eA~ exp —uS^ ■ x it is moving parallel to A"^. The sense of description 
of the ellipse is from the tip of x + eA+exp —ujS~ ■ x towards the tip of 
X -|- eA~ exp —uS~ ■ x that is, from the tip of A+ towards the tip of A~. 

It has been pointed out |1] that the slowness S may not be prescribed a 
priori. Rather, S is written as 

S = NC = Te'^{mm + in), (2.2) 



where 

C = mm + m, m > 1, m ■ fi = 0, |m| = |n| = 1. (2-3) 

In the directional elhpse approach |1], what is prescribed is the directional 
bivector C, or equivalently, the directional ellipse associated with C that is, 
(ill ■ x)^/m^ + (n ■ x)^ = 1, 111 X n-x = 0. On choosing C, then A^, and thus 
S = A^C, are determined from the equations of motion, and A follows as 
an eigenbivector of the acoustical tensor |3]. Let p = m x n be the unit 
normal to the plane of the directional bivector C. There is an infinity of 
choices of directional ellipses, first by varying the magnitude of m, keeping 
111 and h held fixed and orthogonal, and then by rotating (m, n) into another 
orthogonal pair (m*, n*) (say) with p = m* xn* and repeating the procedure. 
Finally p is rotated and the process continued. In this way every possible 
inhomogeneous plane wave solution is obtained. 

For a 'longitudinal' inhomogeneous plane wave the displacement field eu 
may be written 

u= [(5Sexpzw(S-x-t)] + , (2.4) 

where 5 is a constant. Writing 

S=\S\e'^^ S = NC = N{mui + in)=Te'^{mih + ih), (2.5) 

we find 

u = |(5|Texp —uT[{sin (j))'mm ■ x + {cos (f))h ■ x] 

X {mm cos[a;T(cos (f))m,m ■ x — a;T(sin 0)n ■ x — wt + ■?/' + 0] 

— iisin[u;T(cos0)mm ■ x — a;T(sin0)n ■ x — cjt + ■?/' + 0]}. (2.6) 

The planes of constant phase are 

(cos (j))mm ■ X — (sin 0)n ■ x = constant, (2.7) 

and the planes of constant amplitude are 

(sin 0)r7iiii ■ X + (cos 0)n ■ x = constant. (2.8) 

The particle paths are ellipses with principal axes along m and n. We note 
that 

u(0, t) = |(5|T[miiicos('?/' + — ut) — hsm{ip + <P — ^t)], 

u(0, wt = ^ + 0) = |5|Tmiii, 

u(O,wt = ^ + + 7r/2) = |5|Tn. (2.9) 



Hence the sense of description for the motion of the particle initially at 
X = and now on the ellipse fl2.6p with center at x = 0, is from the tip of 
the principal axis along m towards the tip of the principal axis along h. This 
is clearly the case for the motions of all other particles also. 

For circularly polarized waves we take m = 1 so that the planes of con- 
stant phase are orthogonal to the planes of constant amplitude. Also, 

u= |5|Texp — {a;T[(sin(/))m + (cos</))n]-x} 

X [mcos(/t — wt) — nsin(/t — cut)], (2.10) 

where 

K = u;T[(cos0)m — (sin0)n]-x + ip + (p. (2-11) 

The radius of the circle of polarization for the particle initially at x is 

|e(5|Texp — {a;T[(sin0)m + (cos0)n]-x}. (2.12) 

The displacement eu where u is given by O2.10p corresponds to an infinite 
train of circularly polarized longitudinal inhomogeneous plane waves. 



3 Basic equations. 

3.1 Hyper elastic materials. 

We consider a body made of homogeneous isotropic elastic material, initially 
in a state of rest Bq. When the body is deformed to a state B (say), a material 
particle which was at X in Bq moves to the position x, at time t, where 

x = x(X,t). (3.1) 

The left Cauchy- Green strain tensor B is defined by 

B,, = {dx,/dXA){dxj/dXA), (3.2) 

in a rectangular Cartesian coordinate system, fixed in B. 
Three principal invariants of B are I, II, III defined by 

I = tr B, 211 = f - tr (B^), III = det B. (3.3) 

The constitutive equation for the material, relating the Cauchy stress 
tensor T with the strain tensor B is [6] 

T = iVol + NiM - iV_iB"\ (3.4) 



where the material coefficients A^r (r = 0, ±1) are functions of I, II, III. In 
the case of a hyperelastic material, the strain energy density function W, 
measured per unit mass, is characteristic of the material and a function of 
I, II, III: 

1^ = 1^(1,11,111), (3.5) 

and then Nq, Ni, N_i are given by [6j 

1 dW 1 dW 

1 dW 1 dW 

iVi = 2111-2-^, iV_i = 21112^. (3.6) 

3.2 Finite static pure homogeneous deformation. 

We assume that the body is ffist subjected to a finite triaxial static homoge- 
neous deformation, as it is deformed from Bo to B. Then, the corresponding 
three principal axes of deformation constitute a fixed rectangular Cartesian 
coordinate system (0,i,j,k) say, in which the deformation (13. ip is written 
as 

X = AiXi + AaFj + hZk, (3.7) 

where the A's are the principal stretch ratios, assumed, without loss of gen- 
erality, to be ordered as 

Ai > A2 > A3. (3.8) 

In this case, the left Cauchy-Green tensor reduces to 

B = A^i®i + A^j®j + A^k®k, (3.9) 

with corresponding invariants 

l = Xl + Xl + Xl, ll = \l\l + \lXl + XlXl, lll = \l\l\l. (3.10) 

The Cauchy stress necessary to maintain the basic finite static homogeneous 
deformation (13. 7p has components 

To/3 = 0, a 7^ /3 

T„« = iVo + iViA2-iV_iA-2, (3.11) 

where the A^r (r = — 1, 0, 1) are defined as in (13. 6p . at values of I, II, III given 

by dsmD. 



3.3 Superposed 'longitudinal' inhomogeneous plane 
waves. 

Let the body be subjected to a further deformation bringing it from the 
state B to the state B, corresponding to the propagation of an infinitesimal 
'longitudinal' inhomogeneous plane wave of complex exponential type in the 
deformed material, such that the current position x in B, of a particle at x 
given by (13.71) in B, and at X in ^Bq, is given by 

x = x + e{Se*"(^-^-*)}+. (3.12) 

Here, e is a small parameter, such that terms of order higher than e may be 
neglected in comparison with first order terms. 

Now, corresponding to the superposition of the small-amplitude motion 
(I3.12P upon the large deformation (13. 7p are the quantities B, B ,1, II, III, 
A^_i, A^o, -^1, and T, that is, respectively, the left Cauchy-Green tensor 
and its inverse, the first three invariants, the three material constitutive 
parameters, and the Cauchy stress tensor. Each of these quantities may be 
expanded up to the first order in e, and written in the form, 

M = M + e{iujM*e'''^^-''-''> }+,... 

I = tr B = I + e{iwre^'^(^-^-*)}+, . . . 

T = T + e{iujTe'''^^-''-''>}+. (3.13) 



Here, we have. 



B* = S O BS + BS (g) S, 

(B-i)* = -S (g) B-^S - B-^S ® S, 

r = 2(S-BS), 

ir = 2[H(S • S) - IH(S-B-iS)], (3.14) 

Iir = 2IH(S ■ S), 

N* = riVr,i + iriVr,n + IiriVr,in, (F = -1, 0, 1), 



r = N;1 + N;M - iV*iB-^ + A^iB* - iV_i(B-^)*. 
Finally, the mass density p in the current configuration is 

■p = Iir 2p(l - e{iw(S • S)e'"(^-"-*)}+). (3.15) 



3.4 Equations of motion. 

In the dynamic state B, the equations of motion, written in the absence of 
body forces, read 

Up to the first order in e, they are equivalent to T*-S = pS, or 

es + $BS + rB"^s = o, (3.17) 

where 

:= N* + iVi(S-BS) + iV_i(S-l"^S) - p, 

(!>:=N* + Ni{S-S), 

r:=iV!i-iV_i(S-S). (3.18) 

When referred to axes i, j, k, the principal axes of B, Eqs. (13.171) may be 
written 

(6 + <I>A^ + rA~2)C„ = 0, a = 1, 2, 3, no sum, (3.19) 

where C = Cii + C2J + C^k. These equations have to be satisfied for all 
choices of Ci, C2, C3, and for all positive choices of Ai, A2, A3 satisfying 
(13:81) . Noting that 

1 A? Ar^ 



1 A2 A2 
1 A^ Xf 



{Xl- Xl){Xl- Xl){Xl- Xl) /III ^0, (3.20) 

it follows that we must have 



= $ = r = O. (3.21) 

Thus, writing S = A^C, we have 

pN-^ = (2iVo,i + A^i)C-BC + (-2IIIAro,n + A^^OC-B-^C 

+ (2IIiVo,ii + 2IIIiVo,iii)C-C, (3.22) 

and also 

2A^i,i(C-BC) - 2IIIiVi,n(C-B-iC) 

+ [A^i + 2IIAri,n + 2IIIAri,ni] (C ■ C) = 0, (3.23) 

2A^_i,i(C-BC) - 2IIIAr_iji(C-B-^C) 

+ [-iV_i + 2IIAr_i,n + 2lllN_,ju] (C ■ C) = 0. (3.24) 

Here C is assumed prescribed. Then, S = A^C is determined from (13.221) . 
The response functions Ni, N_i must be such that (13.231) and (13.241) are 
satisfied for all choices of C and all positive choices of Ai, A2, A3. 
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4 Propagation of 'longitudinal' inhomogene- 
ous plane waves for any directional bivector 
C. 

Here we find the most general form of the stored energy density for which 
infinitesimal longitudinal inhomogeneous plane waves may propagate in the 
finitely deformed material for any choice of directional bivector C. 

4.1 'Longitudinal' waves of general polarization. 

Equations fl3.23p and f l3.24p must be satisfied for any C. In particular, choose 
C 7^ such that C ■ C = 0, C-B-^C = 0. In that case, C-BC 7^ and it 
follows from ( K23\\ . (Km that 

iVi,i = iV_i,i = 0. (4.1) 

Next choose C 7^ such that C • C = 0, C-IC = 0. In that case, C-B-^C ^ 
and it follows that 

iVi,„ = iV_i,n = 0. (4.2) 

Thus 

A^r = Arr(III), T = -1, 1. (4.3) 

The Nr are functions of III alone, and equations (13.231) and (13.241) reduce 
to 

[A^i + 2IIIA^(] (C ■ C) = 0, [-A^_i + 2IIIiV: J (C ■ C) = 0. (4.4) 

These equations must be satisfied with any C, in particular with C such that 
C-C^O (e.g. C = i + 2y). So, 

A^i + 2IIIA^( = 0, -A^_i + 2IIIA^:i = 0, (4.5) 

and by integration, 

Ni = 2z/Iiri A^_i = 2/iIIli (4.6) 

where ly and /j, are two constant parameters, independent of the strain invari- 
ants. The corresponding hyperelastic material has the following constitutive 
equation, by (13. 4p . 

T = No{l, II, 111)1 + 2z/Iir 2B - 2/iIIl2B~\ (4.7) 



In order that this expression should correspond to that of an hyperelastic 
material, it is seen that A^q has to be independent of I and linear in II, thus 

/ 1 1 d/i \ 

No = 2 ( II 111-2^ + ni2 ^^ J , (4.8) 

so that the strain energy density is given by 

W = ul + i2ll + h{m), (4.9) 

characteristic of the Hadamard material [2j . Here h is an arbitrary function 
of III. We conclude: 

The most general material in which 'longitudinal' inhomogeneous 
plane waves may propagate for any choice of the directional bivec- 
tor C, when it is held in any state of finite static pure homoge- 
neous deformation, is the Hadamard material. 

This same conclusion was reached by John [2\ for the propagation of finite 
amplitude homogeneous plane waves. 

Using (I3.22P and (14.61) . we note that the slowness S corresponding to the 
directional bivector C is given by S = NC, where 

pN-^ = 2{Noj + z/Hr 2)C-1C + 2(-IIIA/'o,ii + /iIH2)C-B"^C 

+ 2(IIAro,„ + IIIAro,„i)C-C. (4.10) 

In the particular case of a hyperelastic Hadamard material, on using (14. 8 p 
and (14. 9p . equation (I4.10p becomes 

pN-^ = 2z/Iir^C-BC - 2/iIH^C-B"^C 



1 1 



+ 4[/iH Iir 2 + IH(IH2/i7]C • C. (4.11) 

For any choice of the directional bivector C, whether it is taken to be lin- 
ear, elliptic, or circular, the corresponding wave train is linearly, elliptically, 
or circularly polarized, respectively, the displacement field being given by 

x = x + e{ArCe'"(^^-"-*)}+. (4.12) 

wherein A^ is given by (l4.1Up . 
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Example 

As an example we may choose C such that 



CI /A^ 


-k) 




CI [xi-xiJ 


fxl 


Then 




C-BC: 


= C-B"^C = 0, 




and (4.111) p;ives 











(4.13) 
(4.14) 



pN'' = -4[/xll lirl + lll(lllk01 ^^' ^'\[f^, ^^l^' ^'^ Cl (4.15) 

^3(^1 + A2) 

Thus NC3 is determined in terms of Aq,, h', and /i, and from (I4.13p . NCi 
and NC2 are also determined. Thus the displacement field 04. 12p may be 
written explicitly in terms of the basic deformation stretches Aq, and fj,, h' . 

We note in passing that (I4.14p may be interpreted |5] in terms of the 
two special central planes for which the sections of the ellipsoids x-Bx = 1, 
x-B^^x = 1, are a pair of similar and similarly situated ellipses. The plane 
of C must coincide with either of the two special central planes and also the 
ellipse of C is similar and similarly situated to the elliptical sections of the 
B and B~^ ellipsoids by the plane of C. 

4.2 'Longitudinal' waves of circular polarization. 

Now we determine the most general material for which infinitesimal 'longitu- 
dinal' inhomogeneous plane waves of circular polarization may propagate for 
any choice of isotropic directional bivector C, when the material is held in 
an arbitrary state of finite static pure homogeneous deformation. Of course, 
as we have seen already, such waves may propagate in a deformed Hadamard 
material. 

The condition that such waves may propagate is that equations 03.23P 
and (I3.24P be satisfied for all isotropic C: C ■ C = and for all positive 
choice of Ai, A2, A3. When C ■ C = 0, equations (I3.23P and (I3.24p reduce to 

^'r,i(C-BC) - llIA/'r,ii(C-B-iC) = 0, {V = -1, 1). (4.16) 

Because Ai, A2, A3 are arbitrary, it follows that 

iVr,i = iVr,n = 0, (r = -l,l), (4.17) 

so that the corresponding constitutive equation is 

T = ^"0(1, II, 111)1 + A^i(III)B - A^_i(III)B"^ (4.18) 
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The corresponding materials are called 'Generalized Hadamard materials'. 
The class encompasses Hadamard materials. 

For any choice of isotropic C, the wave train is circularly polarized, with 
the plane of C being the plane of polarization, and the corresponding slowness 
S = A^C being given through 

pN'^ = (2iVo,i + iVi)C-lC + (iV_i - 2IIIiVo,ii)C-l"^C, (4.19) 

on using (IX^ . 

If the Generalized Hadamard material is to be hyperelastic, then it may 
be shown that Nq must have the form 

A^o = I IH2(IH2 A^i)' - H Iir 2(HI2 Ar_i)' + IH2/i'(III), (4.20) 

and the corresponding form of the strain energy density W is 

2W = I IH2iVi(III) + H HI" 2iV_i(III) + /i(III), (4.21) 

where A^^i, A^_i, and h are arbitrary functions of HI. 

Remark: A universal relation. 

Using equation (14. 19p we derive a universal relation among the wave slow- 
nesses corresponding to three choices of the directional bivector C. We write 
Ci = i + ij, C2 = j + «k, C3 = k + zj. Then Cq, ■ Cq = 0, and 

Oi-]BOi ^ Xi — A2, Oi'IB Oi ^ A]^ — A2 , 

Cs'IBCs ^ A3 — A]^, Cs'IB C3 ^ A3 — -^2 • (4.22) 

Let the value of A^ corresponding to Cq, be N{Ca)- Then 

pN~\Ci) = {2Noj + N,){Xl - Xl) + (Ar_i - 2IIIAro,n)(Ar' - X^^), (4.23) 
etc., and we have the universal relation 

N^^Ci) + iV-2(C2) + N-^iCs) = 0. (4.24) 

5 Propagation of 'longitudinal' inhomogeneous 
plane waves, circularly polarized in special 
planes 

Now we seek elastic isotropic materials which allow the propagation of two 
very special circularly polarized 'longitudinal' inhomogeneous plane waves, 
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irrespective of the basic static homogeneous deformation. The circle of po- 
larization of the waves is to be one or the other of the two central sections of 
the ellipsoid x-B"x = 1 (the B"-ellipsoid), where B is the left Cauchy-Green 
strain tensor fl3.9p corresponding to the basic static deformation. 

We first consider the cases where n = 1 and where n = — 1, and then 
treat the general case. 

5.1 Circle of polarization lies in a plane of central 
circular section of the B-ellipsoid or of the 
B i-ellipsoid. 

The B-ellipsoid is the surface described by x-Bx = 1. Because B is a posi- 
tive definite tensor, there exist two central planes which cut the ellipsoid in 
circular sections. Boulanger and Hayes [5] have proved that the circles are 
described by the bivector A (say) such that A ■ A = A-BA = 0. Hence, the 
wave (I3.12p . with a bivector C satisfying 

C-C = C-BC = 0, (5.1) 

is a 'longitudinal' inhomogeneous plane wave, circularly polarized, with a 
central circular section of the B-ellipsoid as circle of polarization. When 
fl51|) is satisfied, it follows that C-B-^C ^ unless C = 0, so that (K2^ . 
(K2^ reduce to 

iVi,ii = iV-i,ii = 0. (5.2) 

These must hold for all possible positive Ai, A2, A3. Hence A''i = A''i(I,III), 
A^_i = A^_i(I,III) and the corresponding constitutive equation is 

T = No{l, H, 111)1 + A^i(I, III)B - A^_i(I, III)B-^ (5.3) 

This is, of course, inclusive of the Generalized Hadamard materials class 
dug). When (El]) holds, then, from ([S22D, 

pN-^ = (-2niA/'o,ii + iV_i)C-B-iC. (5.4) 

To determine the possible C, we recall that the Hamiltonian decomposition 
of the B tensor is 

B = A^l + |(A? - A^)[h+ ® h- + h- ® h+], (5.5) 

where 

h^ = 6i±4>k, 5 = J|^, 0= j|^, 6' + <p' = l, (5.6) 
6'Xl + 0^A? = Xl 6'Xf + <p'X^' - Xf = 5^\Xl - Xlf /m. (5.7) 
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By dSII]) it follows, using (ES]), that (C ■ h+)(C ■ h") = 0. Hence either 
C ■ h+ = or C ■ h^ = 0, so that the only possible C are 

C = (pi±ii-6k, or C = 01 ± zj + (5k. (5.8) 

For the first pair of possible C, C = 01 + zj — 5k and C = 01 — zj — 5k, 
the polarization circles lie in the same plane (with normal h"^) of central 
circular section of the B-ellipsoid. However the two circles are described in 
opposite senses. A similar comment applies to the second pair of possible C, 
C = 01 ± ij + 5k. The reason why there are four possibilities for C is that if 
C is a solution of (15.11) then so also is its complex conjugate C, because B is 
real and there are just two central circular sections of the B-ellipsoid. 

Now, for all these four possible choices of C, C-B"^C = 5^Aj^ + 0^A]"^ - 
A-2 = 5V^(A? - A^)Vni and so A^ is given by 

pA^-2 _ ^2^2^^2 _ xl)Hir\N^i - 2111^-0,11). (5.9) 

We note that C and A^ are determined completely by the basic deformation. 
The planes of constant phase and the planes of constant amplitude are or- 
thogonal. From (15. 9p . A^ is either purely real or purely imaginary. When A^ 
is real (imaginary), the planes of constant phase (amplitude) are (px ^ 6z = 
const, and the planes of constant amplitude (phase) are y = const. Of course, 
any scalar multiple of S is a possible amplitude bivector. The displacements 
corresponding to (l5.8P i are 

X = X + ee-"^^{(01 ± ij - 5k)e'"[^(^"-^")^*]}+. (5.10) 

Essentially, there are two possible circularly polarized longitudinal inhomo- 
geneous plane waves for which C-BC = 0. 

Also, we note that for (15. 3p to represent a hyperelastic material, the re- 
sponse function A^_i(I, HI) must be independent of I. The response functions 
are then of the form 

N^i = G'(III), A-Q = IIG"(III) + [II/(2III)]G'(III) + K{1, HI), 

A^i = iir 2L(i) + iir 2pir(i, iii)/9i iir 2diH, (5.11) 

with corresponding strain energy density W given by 

2W = jL{l)dl + JK{1, III)Iir 2dIH + H Iir 2G(HI), (5.12) 

where G, K, L are arbitrary functions of their arguments. 
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Similarly, the wave ( I3.12p . with a bivector C satisfying 

C.C = C-B-^C = 0, (5.13) 

is a 'longitudinal' inhomogeneous plane wave, circularly polarized, with a 
central circular section of the B~^-ellipsoid as circle of polarization. When 
fl513|) is satisfied, then fl3:23|l and (KM\ reduce to 

iVi,i = iV_i,i = 0, (5.14) 

which are to be satisfied for all positive Ai, A2, A3. The corresponding con- 
stitutive equation is 

T = Noil, II, 111)1 + A^i(II, III)B - iV_i(II, 111)1"^ (5.15) 

Again, the corresponding class of materials is inclusive of the Generalized 
Hadamard materials class. When (I5.13P holds, then 

piV-2 = (iVi + 2iVo,i)C-BC. (5.16) 

To determine those C for which C ■ C = C-B^^C = 0, we use the Hamil- 
tonian decomposition of the B~^ tensor: 

B-i = \fl + i(Ar' - A^2)[i+ ^ 1- + 1- g) 1+]^ (5.17) 

where 



l± = 5i±0k, 6=J ^_,_^_„ 0=J^A__A_, 6' + <p' = l, 
6^Xf + ^'Xf = \f, 6^X1 + ^^A? - A^ = 6^\Xf - Ar')'lII. (5.1^ 



Here the only possible C are those for which C ■ C = 0, namely C ■ 1^ = 
or C ■ 1^ = 0. Thus, suitable C are 

C = 0i±y-(5k, or C = 01 ± ij + (5k. (5.19) 

For all these four possible choices of C, the corresponding complex scalar 

slowness A^ is given by 

pN~^ = 6^^'^{Xf - Xf)Hll{Ni + 2111^-0,1). (5.20) 

In order that the constitutive equation (I5.15P represent that of a hyper- 
elastic material, the response functions No, Ni, N^i must satisfy certain 
compatibihty equations [3 §86]. Using those compatibility equations it is 
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found that for ( 15.151) the response function A'^i(II,III) must be independent 
of II. Then integrating the compatibility equations it is found that 

iVi = R{111), No = I IIli^[i?(III)IIl2] + S(II,III), (5.21) 

where R and S are arbitrary functions of their arguments and A^_i is a 
solution of 



The characteristics of this equation are 

dlll dll diV_i 



(5.23) 



III II ^N_i + UldS/dir 

the solutions of which depend upon 5*. To make progress we assume 

S = M(III)II + T(III), (5.24) 

where M and T are arbitrary functions of III. In this case the general solution 
of (K23\\ is 

/(II/III, III2 A^_i - /III2 Af (Ill)dlll) = 0, (5.25) 

where / is an arbitrary function, or 

A^_i = -Iir 2/IIl2M(III)dIII + Iir 2/i(II/III), (5.26) 

where h is an arbitrary function. It is found that h must be zero for a 
hyperelastic material. We obtain 

2W = I IIl2i?(III) - II Iiri/IIl2M(III)dIII + /Iir 2T(III)dIII, (5.27) 
and 

A^i = /?(III), A^_i = -Iir 2/IIl2M(III)dIII, 

No = im(III) + I IIIi?'(III) + IIM(III) + T(III), (5.28) 

where R, M, T are arbitrary functions of III. Of course alternative choices 
of 5" will lead to alternative forms for W. 

Finally, we recap for the constitutive model (15. 3 p (or (I5.15P ) that amongst 
the waves which may propagate in the material when it is held in an arbitrary 
state of finite static homogeneous deformation are two infinitesimal circularly 
polarized longitudinal plane waves whose circles of polarization lie in the 
planes of the central circular sections of the B-ellipsoid (or the B~^-ellipsoid) 
corresponding to the finite static deformation. 
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5.2 Circle of polarization lies in a plane of central cir- 
cular section of the B"-ellipsoid. 

Here we determine the most general form of the response functions A^o^ -^15 
A^_i such that the circle of polarization lies in a plane of central circular 
section of the B^-ellipsoid, x-B^x = 1, where B is the left Cauchy-Green 
strain tensor corresponding to the finite static homogeneous deformation. 
As in the previous cases [n = ±1), it is found that the response functions 
A^^i, A^_i depend just upon two invariants of B. They take the forms N-p = 
A^r(III;tr B"). Up to a sign there are four directional bivectors C such that 
two infinitesimal circularly polarized longitudinal plane waves may propagate 
and whose circles of polarization lie in the planes of central circular section 
of the B"-ellipsoid. 
Now we have 

C ■ C = C-B"C = 0, (5.29) 

and equations ([S22D, (IS2SD, (IH^ become 

pA^-2 = (2iVo,i + A^i)C-BC + (-2IIIiVo,ii + N_i)C-n-^C, (5.30) 

and 

Ni^i{C-MC) = IIIA^iji(C-B-^C), 

iV„i,i(C-BC) = IIIiV_i,n(C-B-iC). (5.31) 

Now it may be shown (Appendix) that when 05.291) hold, then 

(C-BC)/(C-B-^C) = III(5trB"/5II)/(9trB"7ai), (5.32) 

so that (15.311) become 

9trB" 9trB" 

^r,^jp - iVr,ii^p = 0, r = ±l. (5.33) 

This suggests that we write A^r(I, II,III) in terms of the set of invariants 
I,III, trB" (n 7^ 1) which is equivalent to the set (I, II, III). We write 

A^r = iVr(I, III, trB"), T = ±1, (5.34) 

and 

dNr dNr dNr <9trB" 
+ 



dl dl 9trB" dl 

dNr dNr dNr <9trB'^ 



dll Oil 5trB" dll 
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(5.35) 



so that (I5.33P becomes 

dNr 



0, r = ±1. (5.36) 



dl 

and thus, because these are to be vahd for all Ai, A2, A3 > 0, we have 

iVr = iVr(III,trB"), T = ±1. (5.37) 

Hence, those materials such that when they are in any state of finite static ho- 
mogeneous deformation, two infinitesimal 'longitudinal' inhomogeneous cir- 
cularly polarized plane waves may propagate, the circle of polarization being 
in the planes of central circular sections of the B"-ellipsoid, where B is the 
strain tensor associated with the finite static homogeneous deformation, have 
constitutive equation 

T = iVo(I, III, trB")l + iVi(III, trB")B - iV_i(III, trB")B-^ (5.38) 

This result is in accord with what we found previously when n = ±1. Indeed, 
forn = 1, trB" = I and fl538D becomes ([53]), whilst forn = -1, trB" = II/III 
and (15.381) may be written in the form of (15.151) . 

Returning to the general case, upon using (15.301) . (15.321) . and (15.351) . we 
find 

pAr-2 = (2—^ + iVi)C-BC + iV_iC-l-^C, (5.39) 

Thus, if C is chosen so that C ■ C = C-B"C = 0, then a circularly polar- 
ized plane wave with slowness S = A^C (A^ given by (I5.39P ) may propagate in 
the material with constitutive equation (I5.38P when it is held in an arbitray 
state of finite static homogeneous deformation. The possible choices of C 
satisfying C ■ C = C-B"C = are 

C = 0„i ± ij - 6nk, C = 0„i ± ij + (5„k, (5.40) 

where 



I \2n \2n / \2" \2n 

*« = (/%3^. *» = y|rr^. « + *S = i. (5.41) 

We note that the terms occuring in (15.391) for all four possible choices of C 
are given by 

C-BC = (A? - Xim - (Xl - Xl)6l 

C-B-^C = [XliXl - Xl)6l - XliXl - Xl)<Pl]/m. (5.42) 
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Special Case n = ^ 

The expressions simplify greatly in the case when n = |, that is, when we 

assume that the circles of polarization lie in the planes of central circular 

1 1 

section of the B2-ellipsoid. Thus C is such that C ■ C = 0, C-B2C = 0, so 

that the possible C are given by 

a/Ai - A3C = a/A2 - Agi ±ij + a/Ai - Aak, 

VAi - A3C = VA2 - Aai ± y - xAT^k. (5.43) 

The corresponding constitutive equation is 

T = iVo(l, 11, 111)1 + iVi (111, trB2)B - iV_i(lll, trB2)B-\ (5.44) 

Using 05. 39 p . the slownesses S = A^C of all four possible circularly polarized 
waves with directional bivectors C given by (I5.43P are such that 

pN-^ = (A, - X,){X, - A3)[2^ + N, + AiA2 + A.A3 + A3A, ^_^j^ ^^^^^^ 



di (A1A2A 



3; 



Appendix 



Here we present a proof of (15.321) . 
For convenience we write 

a = Xl 13 = Xl 7 = A^, (E.46) 

so that 

B = diag (a,/3,7), 1 = a + /3 + 7, 111 = a/37, 

II = a/3 + /37 + 7a, trB" = a" + /3" + 7". (E.47) 

It may be checked that 
da a(/3- 7)111 da -(/3- 7)111 



dl /37(a-/3)(/3-7)(7-«)' dll /3^{a - /3){/3 - ^){^ - a)' 

(E.4J 
and that 



(9trB" _ n[{/3 - 7)a"+i + (7 - a)/3"+i + (a - /3)7"+i] 

(E.49) 



91 (a-/3)(/3-7)(7-a) 

9trB" -n[{f3 - 7)a" + (7 - a)/3" + (a - /3)7"] 



911 (a-/3)(/3-7)(7-a) 
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Thus 


we have the identity 












If C ] 


s such that C ■ C = C-1"C = 0, i.e. 








Cl + Cl + Cl = 0, C^a" + C2/3" + Cl 


7" = 0, 


then 




Cl /3"-7" Cl 7" -a" 
Cl a^'-P''' Cl «"-/?"' 




so that 








C-BC 


Cfa + CIP + Ch 






c-i-ic 


Cl/a + CI//3 + Clh 








jjj (/3 -7)«" + (7 -«)/?" + («- 


/3)7" 




(/3 - 7)a"+i + (7 - a)/3"+i + {a 


- /3)7«+i 


Thus 


when C 


■ C = C-B"C = 0, we have the resuh 
C-BC dtiW/dll 





(E.50) 

(E.51) 
(E.52) 



C-B-iC 9trB"/ai 



(E.53) 



(E.54) 
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